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Evidence of 0(3770) decays to a non-DD final state is observed. A total of 11.8 ± 4.8 ± 1.3 
0(3770) — > J/i{j-k + -k~ events are obtained from a data sample of 27.7 pb~ taken at center-of- 
mass energies around 3.773 GeV using the BES-II detector at the BEPC. The branching fraction is 
determined to be BF(0(377O) -> J/07r + 7r~) = (0.34 ± 0.14 ± 0.09)%, corresponding to the partial 
width of r(^)(3770) -> J/^tt+tt") = (80 ± 33 ± 23) keV. 

PACS numbers: 



I. INTRODUCTION 

The 0(3770) resonance is believed to be a mixture of 
the l 3 Z?i and 2 3 Si states of the cc system pj. Since 
its mass is above the open charm-pair threshold and 
its width is two orders of magnitude larger than that 
of the ip(2S), it is thought to decay almost entirely to 



pure DD |2(. However, Lipkin pointed out that the 
'0(3770) could decay to non— DD final states with a 
large branching fraction There are theoretical cal- 
culations 1^, |3, |6|, 13 that estimate the partial width 
for T (t/>(3770) — > J/-07r + 7r~) based on the multipole ex- 
pansion in QCD. Recently Kuang Q used the Chen- 
Kuang potential model to obtain a partial width for 
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-0(3770) -> J/ipnn in the range from 37 to 170 keV, 
corresponding to 25 to 113 keV for -0(3770) — > J/07r + 7r _ 
from isospin symmetry. In this paper, we report evidence 
for 0(3770) — > J/ipir + iT~~ based on a data sample of 27.7 
pb _1 taken in the center-of-mass (cm.) energy region 
from 3.738 GeV to 3.885 GeV using the upgraded Beijing 
spectrometer (BES-II) at the Beijing Electron Positron 
Collider (BEPC). 

II. THE BES-II DETECTOR 

The BES-II is a conventional cylindrical magnetic de- 
tector that is described in detail in Ref. [8|. A 12-layer 
Vertex Chamber (VC) surrounding the beryllium beam 
pipe provides input to the event trigger, as well as co- 
ordinate information. A forty-layer main drift chamber 
(MDC) located just outside the VC yields precise mea- 
surements of charged particle trajectories with a solid an- 
gle coverage of 85% of 4-zr; it also provides ionization en- 
ergy loss (dE/dx) measurements which are used for parti- 
cle identification. Momentum resolution of 1.7%yl + P 2 
(p in GeV/c) and dE/dx resolution of 8.5% for Bhabha 
scattering electrons are obtained for the data taken at 
■s/s = 3.773 GeV. An array of 48 scintillation counters 
surrounding the MDC measures the time of flight (TOF) 
of charged particles with a resolution of about 180 ps for 
electrons. Outside the TOF, a 12 radiation length, lead- 
gas barrel shower counter (BSC), operating in limited 
streamer mode, measures the energies of electrons and 
photons over 80% of the total solid angle with an energy 
resolution of a E /E = 0.22/ ' \[~E (E in GeV) and spatial 
resolutions of = 7.9 mrad and az = 2.3 cm for elec- 
trons. A solenoidal magnet outside the BSC provides a 
0.4 T magnetic field in the central tracking region of the 
detector. Three double-layer muon counters instrument 
the magnet flux return and serve to identify muons with 
momentum greater than 500 MeV/c. They cover 68% of 
the total solid angle. 

III. DATA ANALYSIS 

A. Monte Carlo simulation 

To understand the main source of background in the 
study of the decay 0(3770) — > J/07r + 7r~, we developed 
a Monte Carlo generator. The Monte Carlo simulation 
includes the initial state radiation (ISR) at one loop or- 
der, in which the actual center-of-mass energies after 
ISR are generated according to Ref. The 0(25) 

and 0(3770) are generated using energy dependent Breit- 
Wigner functions according to Eq. (38.53) of Ref. [UJ in 
which the ratio of T e i(s)/T to t(s) = F^,(2S)^e+e- ft tot 
and the branching fraction of 0(25) — ► J/07r + 7r~ in 
the formula are assumed to be constant. The beam en- 
ergy spread (<JE bcam — 1-37 MeV) is taken into account 
in the simulation. Since there is no unique description 



and solution for the low energy tttt production ampli- 
tude the correction of the decay rate due to the 
7T7r production amplitude is neglected in the making of 
the event generator. However the effect of the varia- 
tions in the correction to the decay rate on the estimated 
number of 0(25) — > J/07r + 7r _ is considered in the final 
background subtraction in subsection E. Fig. ^ shows 
the distribution of J/07r + 7r~ events with J/0 — > l + l~ 
(I = e or fj) as a function of the actual energy remain- 
ing after initial state photon radiation, which is deter- 
mined by our Monte Carlo generator, where the branch- 
ing fraction for 0(3770) — > J/07r + 7r~ is set to be 0.35%, 
while the branching fractions for 0(25) — * J/07r + 7r~ and 
J/0 — > l + l~ are taken from the Particle Data Group 
(PDG) [13. 

There are two peaks which are around 3.686 GeV and 
3.773 GeV in Fig. 2] where the J/07r + 7r~ events from the 
0(25) decay are given by the dotted histogram, while the 
events from the 0(3770) decay are given by the solid his- 
togram. There are two components in the higher mass 
peak. One is from 0(3770) production, and the another is 
from 0(25) production which is due to the tail of 0(25) 
Breit-Wigner function. This type of 0(25) production 
(called type B of 0(25) in this Letter) is indicated by 
the dotted histogram around 3.773 GeV. The J/07r+7r~ 
events in the lower mass peak are produced around the 
peak of the 0(25) Breit-Wigner function (called type A 
of 0(25) in this Letter), and are due to ISR energy re- 
turn to the 0(25) peak. The "type B" of 0(25) is the 
main source of background events in the experimental 
study of the decay 0(3770) — > J/07r + 7r _ . This back- 
ground event has the same topology as that as the de- 
cay of 0(3770) — > J/07r + 7r~. To get the number of 
0(3770) — > J/07r + 7r~ signal events, the number of the 
background events of 0(25) — > J/07r + 7r _ has to be sub- 
tracted from the observed candidate events of J/07r + 7r~ 
based on analyzing the Monte Carlo events. 

In the energy region from 3.738 to 3.885 GeV in Fig.^ 
there are 1724 0(3770) -> J/07T+7T" events and 747 
0(25) — > J/07r + 7r~ events. The generated events as 
shown in Fig. ^ are put through the full detector sim- 
ulation based on the GEANT simulation package. The 
fully simulated events are used for study of the main 
background. 



B. Events selection 

To search for the decay of 0(3770) — > J/07r + 7r~, 
J/0 — > e + e~ or /i + /i~7r + 7r~ and e + e~7r + 7r~ can- 

didate events are selected. These are required to have 
four charged tracks with zero total charge. Each track is 
required to have a good helix fit, to be consistent with 
originating from the primary event vertex, and to satisfy 
| cos 6*| < 0.85, where is the polar angle. 

Pions and leptons must satisfy particle identification 
requirements. For pions, the combined confidence level 
(CL), calculated for the 7r hypothesis using the dE/dx 
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FIG. 1: The numbers of -0(3770) -> J/^tt+tt" (solid line) and 
ip(2S) — > J/if)-K + ir~ (dotted line) as a function of the actual 
energy remaining after ISR, where J/0> is set to decay to 
l + l~ ; the events are generated with the Monte Carlo generator 
at the cm. energies at which the data were collected from 
3.738 to 3.885 GeV. The insert on the right-top shows the 
distribution of the energy at which tp(2S) — > J/0>7r + 7r~ events 
are produced. 



and TOF measurements, is required to be greater than 
0.1%. In order to reduce 7 conversion background, in 
which the e + and e~ from a converted 7 are misidenti- 
fied as 7r + and n~ , an opening angle cut, 0^+^- > 20°, is 
imposed. For electron identification, the combined con- 
fidence level, calculated for the e hypothesis using the 
dE/dx, TOF and BSC measurements, is required to be 
greater than 1%, and the ratio CL e /(CL e + CL fl + CL 7r + 
CLk) is required to be greater than 0.7. If a charged 
track hits the muon counter, and the z and r0 positions 
of the hit match with the extrapolated positions of the 
reconstructed MDC track, the charged track is identified 
as a muon. 

The candidate events of e + e~7r + 7r~ or /i + /.i~ir + ir~ sat- 
isfying the above selection criteria are further analyzed 
by using two different analysis methods to be discussed 
in subsection C and subsection D. 



C. Analysis of the tt + -k recoil mass 

The mass recoiling against the 7r + 7r _ system is calcu- 
lated using 



Aft,Ec(7r + 7T 



where E cm is the cm. energy, E 7r + 7T - and P^+ v - are 
the total energy and momentum of the 7r + 7r _ system, 
respectively. 

Fig. [3 shows the distribution of the masses recoiling 
against the 7r + 7r~ system for candidate events with total 
energy within ±2.5<te + _ i+i „ of the nominal cm. en- 
ergy at which the events were obtained and with a dilep- 
ton invariant mass within ±150 MeV of the J/tp mass, 
where <te + + is the standard deviation of the distri- 
bution of the energy of the tt + tt~1 + 1~ . Two peaks are 
observed. The higher one is from the "type A" of ip(2S) 
events produced by radiative return to the peak of the 
ip(2S), while the small enhancement around 3.1 GeV is 
mostly from -0(3770) decays, but also contains the con- 
tamination of the "type B" of ip(2S) decays. This is 
confirmed by analyzing the Monte Carlo sample gener- 
ated with the Monte Carlo generator as mentioned be- 
fore. Fig. |3| shows the same distribution from analysis 
of the Monte Carlo events for i>(2S) and 0>(377O) pro- 
duction and decays to J/-07r + 7r~ with J/tp — > final 
states. The Monte Carlo events are generated at the cm. 
energies at which the data were collected. The size of the 
Monte Carlo sample is twenty times larger than the data. 
There are also two peaks; the higher one is from the "type 
A" of ip{2S) events; the small enhancement around 3.1 
GeV consist of two components. One is from the "type 
B" of ip{2S) event production and decays as shown by 
the solid line histogram; the other one is from -0(3770) 
production and decays which come from the events as 
shown by the solid histogram in Fig^ The error bars in 
Fig. |3| are based on the total number of observed events 
of 0^(3770) -> J/-07T+7T- and i/;(2S) -> J/ipir+n-. The 
differences between the error bars and the histogram as 
shown around 3.1 GeV correspond to 0(3770) production 
and decay to J/-07r + 7r~. 

Monte Carlo studies show that the distributions of the 
masses recoiling from the ir + ir~ system for the "type A" 
of 0'(25 l ) events can be described by two Gaussian func- 
tions. One Gaussian function is for the "type A" of ip(2S) 
production from the events for which the nominal cm. 
energy is set at 3.773 GeV, and the other one is from 
the events for which the nominal cm. energies are set 
off 3.773 GeV. Using triple Gaussian functions, one of 
which describes the peak near the J/-0 mass and two of 
which represent the second and the third peaks of the 
events from the "type A" of ip(2S), and a first order 
polynomial to represent the background to fit the mass 
distributions as shown by the solid histograms for both 
the data (Fig. [2J and the Monte Carlo (Fig. [3J sample, 
we obtain a total of 25.5 ± 5.9 J/ip — > signal events 
from both the ^(3770) -> J/^tt+tt" and the "type B" of 
ip(2S) -> J/ipn+n- events and 220.5 ± 26.0 "type B" of 
ip{2S) — > J/077r + 7r~ events, respectively. The curves give 
the best fits to the data and the Monte Carlo sample. 
The fitted peak positions and standard deviations of the 
Gaussian functions used for the fits in Fig. [21 and Fig. |3| 
arc listed in Table I. 
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TABLE I: Summary of the fitted results of the data and the 
Monte Carlo sample in Fig. |2] and Fig. [21 



Figure Peak Mass [MeV] a M [MeV] 

Peakl 3097.8 ± 3.0 9.9 ± 2.4 
H Peak2 3182.4 ± 1.8 23.3 ± 2.8 
Peak3 3182.5 ± 0.6 7.6 ± 0.7 

Peakl 3099.4 ±1.3 13.1 ±2.9 
El Peak2 3180.8 ± 0.4 22.4 ± 0.7 

Peak3 3185.0 ± 0.2 7.4 ± 0.2 
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FIG. 3: The distribution of the masses recoiling from tt + tt~ 
for the Monte Carlo events of il>(2S) — > J/ipn + n~~ and 
1^(3770) — » J/i{j-k + -k~ with J/ip — > £ + Z~; these events are gen- 
erated with the Monte Carlo generator; where the error bars 
represent the sum total of the two components, while the his- 
togram is for the ip(2S) — > J/ip% + -K~ from both the "type A" 
and "type B" of ip(2S); the curves give the best fits to the re- 
coil mass spectrum from the Monte Carlo ip(2S) —* J/ipn + 7v~ 
events only; the insert on the left-top shows the mass distri- 
bution in a local region around 3.1 GeV; here two compo- 
nents from the 0(3770) -> J/ipTT+%' and the "type B" of 
ip(2S) — » J/i{j-k + -k~ are clearly seen. 



FIG. 2: The distribution of the masses recoiling from the 
7r + 7r - system for l + l~iv + iv~ events; the insert on the left-top 
shows the mass distribution in a local region around 3.1 GeV; 
the curves give the best fit to the recoil mass spectrum, see 
text. 



D. Kinematic At 

In order to reduce background and improve mo- 
mentum resolution, candidate events are subjected to 
four-constraint kinematic fits to either the e + e~ — > 
/i + /i~7r + 7r~ or the e + e~ — * e + e~7r + 7r~ hypothesis. 
Events with a confidence level greater than 1% are ac- 
cepted. Fig. ^] shows the dilepton masses determined 
from the fitted lepton momenta of the accepted events. 
There are clearly two peaks. The lower mass peak is 
mostly due to -0(3770) — > J/ipir + n~, while the higher 
one is due to the "type A" of ip{2S) — ► J/07r + 7r~. Since 
the higher mass peak is produced by the radiative re- 
turn to the ip(2S) peak, its energy will be approximately 
3.686 GeV, while the cm. energy is set to the nominal 



energy in the kinematic fitting. Therefore, the dilepton 
masses calculated based on the fitted lepton momenta 
from i/j(2S) — > J/-07r + 7r~, J/ip — > l + l~ are shifted up- 
ward to about 3.18 GeV. 

A maximum likelihood fit to the mass distribution in 
Fig. using three Gaussian functions to describe the 
mass distribution of the 1 + 1~ti + it + combinations and a 
first order polynomial to represent the broad background 
as used to fit the 7r + 7r~ recoil mass distributions in Fig. [21 
and Fig. yields a J/ip mass value of 3097.8 ± 3.0 MeV 
and a signal of 17.8 ±4.8 J/tp —> l + l~ events. The curves 
give the best fit to the data. 

As discussed in Section C, there is a contribution from 
the "type B" of ip(2S) — > J/ipir + ir~ that can pass the 
event selection criteria and can lead to an accumulation 
of the recoil masses of the tt + tt~ and/or the fitted dilep- 
ton masses around 3.097 GeV. This is the main source of 
background to -0(3770) — > J/i(m + Tr~. Fig. shows the 
distribution of the fitted dilepton masses of the Monte 
Carlo events of ^(3770) -> J/^tt+tt" and 0(25) -> 
J/07r + 7r~ with J/0 — > as shown in Fig. ^ Here 

the histogram shows the dilepton mass distribution for 
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FIG. 4: The distribution of the fitted dilepton masses for the 
events of l + l~ir + iv~ from the data; the hatched histogram 
is for fi + n~ir + iv~ , while the open one is for e + e~7r + 7r~; the 
curves give the best fit to the data. 



ip(2S) — ► J/ipTT + n~ only. The higher mass peak is due to 
the "type A" of tp(2S) — ► J/ipTr + ir~ events, and the lower 
one is from the "type B" of tp(2S) — > J/ipTr + ir~ events. 
The error bars show the sum total of the observed events 
of ip(2S) -» J/^tt+tt- and ^(3770) -» J/ipn+n-. The 
differences between the error bars and the histogram cor- 
respond to the observed events of i/>(3770) — ► J/ipn + n~ . 
Fitting the mass distribution for the tp(2S) events only 
(histogram) with the same triple Gaussian functions as 
mentioned before yields 119 ± 12.1 J/ip events from the 
"type B" of ip(2S) -» J/^tt+tt- decays. 



E. Other background and background subtraction 

1. Other background 



FIG. 5: The distribution of the fitted dilepton masses for 
the events of £ + Z~7r + 7r~ from the Monte Carlo sample of 
^(3770) — * J/ip-K + -K~ and ip(2S) — > J/ipn + Tr~ which are gen- 
erated with the Monte Carlo generator (see section A); the 
histogram is for ip(2S) — > J/ipir + -K~ , while the error bars are 
the sum of V(3770) -> J/ipTi+Tv- and ip(2S) -> J/^Tr+Tr^, 
where the J/tp is set to decay to l + l~ . 
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Some physics processes, such as two-photon events, 
e+e~ — > e+e~^+/i~ (where the slow muons are misiden- 
tified as pions) and e + e~ — » e + e~7r + 7r~, e + e~ — > t + t~ 
and e+e~ — > Dl) could be sources of background. 

To check if there are some background contamina- 
tions in the observed J/ipir + ir~ events due to the pos- 
sible sources of background, we generated 1 x 10 5 two- 
photon Monte Carlo events (which is about 4 times larger 
than the data), 6 x 10 5 e + e~ — ► hadrons Monte Carlo 
events (which is about 1.6 times larger than the data), 
and 2.3 x 10 6 e+e~ -> DD Monte Carlo events (which 
is about 13 times larger than the data), where the D 



JUL 



3.1 



3.2 



3.3 



Fitted mass (GeV/c ) 



FIG. 6: The distributions of the dilepton masses of the events 
of l+r-K+K' from the data taken at 4.03 GeV with the BES-I 
at the BEPC. 
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FIG. 7: The distribution of the masses recoiling against 
the tt + -k~ system calculated using the measured momenta 
for events that pass the kinematic fit requirement, where the 
hatched histogram is for events of fi + /i~ir + iv~ and the open 
one is for e + e~ir + -K~ . 

and D mesons are set to decay to all possible final states 
according to the decay modes and branching fractions 
quoted from PDG 10] . These Monte Carlo events are 
fully simulated with the GEANT-based simulation pack- 
age. None of the simulated possible background events 
were misidentified as J/ipTT + ir~ events. 

The candidate J/ipn + Tr~ events could also be produced 
in the continuum process, such as e + e~ — ► l + l~ir + ir~ 
and e + e~ — > t + t _ , and satisfy the selection criteria. 
From analyzing a sample of 6.6 pb _1 taken at 3.65 GeV 
with the BES-II detector, a sample of 5.1 pb _1 taken in 
the energy region from 3.544 to 3.600 GeV and a sam- 
ple of 22.3 pb" 1 taken at 4.03 GeV with the BES-I de- 
tector, no significant J/ipir + ir~, J/tp — ► events are 
observed. Fig.||)]shows the distribution of the fitted dilep- 
ton masses of the events of l + l~ir + ir~ which satisfy the 
selection criteria; these events are from the data taken 
with the BES-I detector at 4.03 GeV. The distribution 
of the fitted dilepton masses is flat, which is consistent 
with the background distribution. Hence the continuum 
background is negligible. 



2. Number of background 

After normalizing to the total luminosity of the data 
set, we estimate that there are 11.0 ± 1.3 ± 2.4 back- 
ground events from the "type B" of ip(2S) — > J/tpir + ir~ 



in the 25.5 ± 5.9 J/ip — > signal events obtained by 
fitting to the it + it~ recoil mass distribution of Fig.[5|and 
6.0 ± 0.5 ± 1.3 background events from the "type B" of 
ip(2S) -> J/ipTr+n- in the 17.8 ±4.8 J/0 -> l+l~ signal 
events obtained by fitting to the fitted dilepton mass dis- 
tribution of Fig. 21 where the first errors are statistical 
and the second are systematic. The later arise from the 
uncertainties (±1.1) and (±0.6) in the ip(2S) resonance 
parameters and the uncertainties (±2.2) and (±1.2) com- 
ing from the ambiguities of the knowledge of the low en- 
ergy 7T7T production amplitude in tp(2S) — > J / ij>iT + tt~ . 
The two terms correspond to the uncertainty on the pro- 
duction amplitute predicted by different theoretical mod- 
els |l l| for analyzing the it + tt~ recoil mass spectrum and 
the fitted dilepton mass spectrum, respectively. The the- 
oretical models are based on the PCAC and current al- 
gebra or chiral pertubative theory predictions in which 
the various 7r + 7r~ rescatering corrections are taken into 
account to get better unitarity behavior at higher energy. 



F. Number of signal events 0(3770) — > J/07T 7T 

The probability that the 17.8 events observed are due 
to a fluctuation of the 6.0 ± 0.5 ± 1.3 events is 1.1 x 
10 -3 . After subtracting the numbers of the background 
events, 14.5 ±6.4 ±2.4 and 11.8 ±4.8 ±1.3 signal events 
of ■0(3770) — * J/ifjTr + 7T~ are retained from analyzing the 
recoil masses of ir + ir~ (see Fig. 2) and the fitted dilepton 
masses (see Fig.0J of the events l + l~ir + ir~ , respectively. 

In this analysis, the possible interference between the 
ip(2S) background and the ^(3770) signal is neglected 
since the decay wave functions of 0(3770) and ip(2S) are 
orthogonal 6]. The BES detector is symmetric enough 
in the spatial direction and there is no bias for the event 
selections about the momentum direction of the particles. 
Therefore the interference terms cancel after integrating 
over the pion momenta. 

To test whether there is any bias in the kinematic fit, 
we examine the tt + tt~ recoil mass distribution for the 
events passing the kinematic fit requirements. Fig. [7| 
shows the recoil mass distribution, where the recoil 
masses are calculated as mentioned in Section C, but 
the events are not required to satisfy the total energy 
cut and dilepton invariant mass cut. There are also two 
peaks, similar to those in Fig. 0] observed clearly. Fit- 
ting to the mass spectrum with the same functions as 
described above yields a J/ip mass value of 3100.1 ± 3.8 
MeV and a signal of 17.2 ±5.0 events, consistent with the 
17.8±4.8 signal events obtained by fitting to the dilepton 
mass distribution in Fig. 

Tabic II summarizes the fitted peak positions and stan- 
dard deviations of the Gaussian functions used for the fits 
in Fig.U Fig. Eland Fig. [7] 



7 



TABLE II: Summary of the fitted results of the data and 
Monte Carlo sample in Fig. 21 Fig.^and Fig. Q 



Figure Peak 


Mass [MeV] 


om [MeV] 




Peakl 


3097.8 ± 3.0 


9.9 ±2.4 


Fig 


. [2 Peak2 


3173.1 ±5.5 


24.8 ±5.7 




Peak3 


3180.9 ±2.4 


9.3 (fixed) 




Peakl 


3098.8 ±0.7 


9.1 ±0.5 


Fig 


. Peak2 


3169.2 ±0.7 


22.2 ±0.5 




Peak3 


3181.8 ±0.2 


9.3 (fixed) 




Peakl 


3100.1 ±3.8 


13.2 ±3.5 


Fig 


. □ Peak2 


3177.1 ±3.9 


19.6 (fixed) 




Peak3 3185.4 (fixed) 


7.0 (fixed) 



G. The number of 0(3770) produced 

The total number of "0(3770) events is obtained from 
our measured luminosities at each cm. energy and from 
calculated cross sections for 0(3770) production at these 
energies. The Born level cross section at energy E is 
given by 

b /pn 127rr ee r tot (£') 

^(3770) \E) - {E 2 „ M 2 )2 + M 2 T 2 ot{E y 

where the ■0(3770) resonance parameters, T ee and M, 
are taken from the PDG 01 an d F tot (E) is chosen to 
be energy dependent and normalized to the total width 
r to t at the peak of the resonance ^(j^^El- ^ n or ~ 
der to obtain the observed cross section, it is neces- 
sary to correct for ISR. The observed 0(3770) cross sec- 
tion, a^ 377Q) (s nom ), is reduced by a factor g(s nom ) = 

^^3770) (^ nom )/ f7 '^(377o) (^ nom ) 3 where s n om is the cm. 
energy squared and 0^(3770) {s nom ) is the Born cross sec- 
tion. The ISR correction for 0(3770) production is cal- 
culated using a Breit-Wigner function and the radia- 
tive photon energy spectrum With the calcu- 
lated cross sections for 0(3770) production at each energy 
point around 3.773 GeV and the corresponding luminosi- 
ties, the total number of "0(3770) events in the data sam- 
ple is determined to be A^/377 = (1.85 ± 0.37) x 10 5 , 
where the error is mainly due to the uncertainty in the 
observed cross section for 0(3770) production. 



IV. RESULT 
A. Monte Carlo efficiency 

The efficiencies for reconstruction of the events 
of 0(3770) — > J/07r + 7r~ with J/0 — > e + e~ and 
J/0 — + /Lt + M _ are estimated by Monte Carlo simula- 
tion. Monte Carlo study shows that the efficiencies 
are e^ 3770 ^j/^+^- t j/^ e+e - = 0.146 ± 0.003 and 



eV(377o)^,//-07r+7 r -„//^-^ AI + Al - = 0.174 ± 0.003, where the 
errors are statistical. These give the averaged efficiency 
for detection of J/0 — > e + e~ and J/0 — > events 
to be e</,(3770)^J/^+ 7 r-,j/^-»;+;- = 0.160 ± 0.002. 

B. Branching fraction and partial width 

Using these numbers and the known branching frac- 
tions for J/0 — * e + e~ and /J + /z~ the branching 
fraction for the non—DD decay "0(3770) — > J/07r + 7r~ is 
measured to be 

55(0(3770) -► J/07T+7O = (0.338 ± 0.143 ± 0.086)%, 

where the first error is statistical and the second system- 
atic arising from the uncertainties in the total number of 
0(3770) produced (20% ), tracking efficiency (2.0% per 
track), particle identification (2.2%), background shape 
(6%), 0(25) -> J/ipn+ir- background subtraction (12%) 
and the averaged branching fraction for J/0 — > 
(1.2%). Adding these uncertainties in quadrature yields 
the total systematic error of 25.5%. 

Using r to t from the PDG ^(j> this branching fraction 
corresponds to a partial width of 

T(0(377O) -> J/07T+7T-) = (80 ± 33 ± 23) keV, 

where the first error is statistical and the second system- 
atic. The systematic uncertainty in the measured partial 
width arises from the systematic uncertainty in the mea- 
sured branching fraction (25.5%) and the uncertainty in 
the total width of 0(3770) (11.5%) [ijj. 

As a consistency check, we can use the number of signal 
events 0(3770) — > J/07r + 7r~ obtained from analyzing the 
recoil masses of the tt + tt~ from the events 1 + 1~tt + it~ to 
calculate the branching fraction. The detection efficiency 
is £v(3770)-»,//^7r+7r-,j/^-*i+i- = 0.194. These numbers 
yield a branching fraction 55(0(3770) — > J/07r + 7r~) = 
(0.342±0.142±0.097)%, which is in good agreement with 
the value obtained above, indicating that the kinematic 
fit result is reliable. 



V. SUMMARY 

In summary, the branching fraction for 0(3770) — > 
J/07r + 7r - has been measured. From a total of (1.85 ± 
0.37) x 10 5 0(3770) events, 11.8 ± 4.8 ± 1.3 non-55 de- 
cays of 0(3770) — > J/07T+7T - events are observed, leading 
to a branching fraction of 55(0(3770) — > J/07r + 7r~) = 
(0.34 ± 0.14 ± 0.09)%, and a partial width T(0(377O) -> 
J/07r+7r-) = (80 ± 33 ± 23) keV. 
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